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The analysis of correlation function data obtained by Monte Carlo simulations of the two-dimensional 
4-state Potts model, XY model, and self-dual disordered Ising model at criticality are presented. We study 
the logarithmic corrections to the algebraic decay exhibited in these models. A conformal mapping is used 
to relate the finite-geometry information to that of the infinite plane. Extraction of the leading singularity 
is altered by the expected logarithmic corrections, and we show numerically that both leading and correction 
terms are mutually consistent. 
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A second order phase transition occurs at very spe- 
cial points in the parameter space of a model, i.e. at a 
fixed point of the renormalization equations. The eigen- 
values of the linearized renormalization equations define 
the scaling dimensions of the corresponding directions 
related to the scaling fields. Positive eigenvalues are 
associated to relevant scaling fields while negative ones 
correspond to irrelevant fields. In some cases, there may 
exist a line of fixed points along which critical expo- 
nents are varying. This occurs when a marginal field 
(with vanishing scaling dimension) is identified in the 
model. The scaling dimensions completely characterize 
the critical properties of the model, e.g. the power laws 
of the physical quantities. As an example, the corre- 
lation functions exhibit an algebraic decay at the fixed 
point, e.g. for a scaling field density </>(r), 



only), discrete spin variables <j w = 1, 2 . . . q are located 
at the sites w of a square lattice (/i specifies the unit 
lattice vector in the two directions) and interact with 
nearest neighbours, 
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'in two dimensions (rj^ should be replaced by d — 2 + 
in arbitrary dimension). In some special cases, this sim- 
ple behaviour is modified by multiplicative logarithmic 
terms, 



The value of the number of states per spin, q, can be 
seen as a parameter of the model. The special case 
q — > 1 corresponds to the percolation problem and ex- 
hibits a second order phase transition. The value q = 2 
corresponds to the Ising model, so when q varies the 
universality class changes. The transition is of second 
order as long as q < 4 and the value q = 4 coincides 
with the end of the line of fixed points. Along this line, 
the exponent of the order parameter correlation function 
varies according to 
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This situation occurs at the end of a line of fixed points 
or when the system is perturbed by a marginally irrel- 
evant operator. Examples in two dimensions are given 
by 4-state Potts model, XY model, or disordered Ising 
model. 

In the case of the g-state Potts model (from now on, 
we consider the case of the two dimensional problem 



where the value of m = tt/cos^ 1 ^^) — 1 parametrizes 
the number of states. At q = 4, the correlation func- 
tion exponent takes the value ri a {q = 4) = j. Above 
q = 4 the transition becomes of first order. The transi- 
tion temperature of the 2d Potts model is exactly known 
from duality requirements [2], K c = ln(l + -Jq). 

The two-dimensional classical XY model is another 
example of such a scenario, although the transition is in 
its very nature quite different. The Hamiltonian now de- 
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scribes the interaction of classical two dimensional unit 
spins, er w = (cos 9 w , sin 6 w ), 

Hxy 



(5) 
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At low temperature, as a consequence of the Mermin- 
Wagner-Hohenberg theorem [3] , there is no finite macro- 
scopic magnetization, but the spin-spin correlation func- 
tion decays algebraically in the so-called quasi-long- 
range-ordered phase. The temperature plays the role 
of a marginal field and the decay exponent T] a contin- 
uously increases up to a limiting value rja (XY) = j 
at a transition temperature named after Berezinskii- 
Kosterlitz-Thouless (BKT) 4 . In this low tempera- 
ture critical phase, ordering is prevented by collective 
excitations, spin waves, and localized excitations, vor- 
tices, which appear in increasing number as the temper- 
ature is increased. These latter topological defects are 
bounded in pairs and the transition to a completly dis- 
ordered phase with exponential decay of the correlations 
is reached at the BKT temperature when unbinding of 
the pairs occurs [S]. The BKT transition thus corre- 
sponds to the end of a continuous line of fixed points in 
the low temperature phase. The transition temperature 
K c ~ 0.893, as well as the expression of the temperature 
dependence of the exponent -q^ are not known exactly. 
The mechanism of the transition may also be under- 
stood from the role of the vortex chemical potential. 
It is a relevant variable in the low temperature phase, 
which becomes marginally irrelevant at the BKT tran- 
sition, thus producing the essential singularities and the 
logarithmic corrections. 

A third example is provided by the random bond 
Ising model in two dimensions. According to Harris 
criterion [Jj], quenched disorder is a relevant variable 
when the specific heat exponent of the pure system un- 
der consideration is positive. In this case, a new fixed 
point leading to a new universality class is expected. In 
the case of the 2d Ising model, since a = in the pure 
model, randomness is only a marginal variable which 
could either produce a continuous variation of the ex- 
ponents with the amplitude of disorder or logarithmic 
corrections to the unchanged leading critical behaviour 
if the disorder is eventually marginally irrelevant. After 
an interesting debate in the eighties |7j, the second sce- 
nario has been recognized to be correct |Hj and the corre- 
lation function exponent keeps its value r\ a (RBIM) = j . 
The Hamiltonian is the one given in equation (0, with 
q = 2 and with random nearest neighbour interactions 
K Wt/x along the bonds between sites w and w + fi. When 
these couplings are taken from a binary probability dis- 
tribution with equal probabilities for both strengths, 



HK W J = \ U W JS{K W ^ - K x ) + 5(K W ^ - K 2 )], the 
critical temperature follows from duality, (<sxp(Kf) — 
l)(exp(i^) - 1) = q (= 2 here). 

All these models have the common property that the 
order parameter correlation function G a is expected, 
from renormalization group arguments, to behave ac- 
cording to equation J5J) with 9 a ~ ^ (XY model) and 
-| (Potts model and RBIM) jHUHIIinj. One may then 
define a local effective exponent r] e g(r) (r = |ri — r2 1), 
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To make this formula more explicit, let us put some 
typical numbers. Suppose that we want to produce nu- 
merical simulations, in order to check equation |2"| . with 
let say a relative accuracy of Arj/rjc = 
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Tj„ lnr 
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The exponent \9 a \ might be estimated to be of the or- 
der of the value of r] a , so that one has to reach values 
of r as large as r ~ exp^A^/^) -1 ) ~ 10 43 ! Even 
with a value of \6 a \ — i, one still needs a sample much 
larger than L > 10 21 in order to reach distances between 
spins which are large enough to become sensitive to the 
presence of the log in equation (J5J). This is definitely 
not possible and that might be the reason why it was 
almost impossible to produce reliable data ^ to corrob- 
orate equation J5J. The strategy is thus to accept the 
leading behaviour (the value of r] a as predicted by RG), 
and to extract the value of 9 a . This is also a difficult 
route, and contradictory results were reported in the lit- 
erature. In the case of the XY model for example, fixing 
rj a = 4 and using Monte Carlo simulations (susceptibil- 
ity and correlation length data), Kenna and Irving 1 1 1 j 
reported a value close to B a ~ 0.046(20), Janke [T2"] 
obtained 9 a ~ 0.054(2) from data at criticality, or 
9 ^ ~ —0.112(4) using high temperature data. Also with 
high temperature data, Patrascioiu and Seiler ob- 
tained 9 a ~ —0.154(92) and from strong coupling ex- 
pansion Campostrini et al. ^3 had results depending on 
the lattice symmetry, 6 a ~ -0.090(6) to -0.084(12) 2 \ 
There are even more controversial results reported by 
Balog et al. ^H] who suggest 9 a = plus additive cor- 
rections, or Kim |16j who is in favor of ordinary scaling 
rather than essential singularities. We are not aware 
of any direct verification of the presence of logarithmic 
terms directly in the correlation function in the case 
of 4-state Potts model or disordered Ising model, but 
logarithmic terms have been found to be compatible 
with finite-size scaling data of susceptibility and specific 
heat ^7] in the 4-state Potts model and in the random 
bond Ising model |18| . 

1 'Both by simulations and by series expansions. 
2 'ln the literature, B a is often referred to as —It. 
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There are conceptual difficulties with the above- 
mentioned approach applied to numerical simulations. 
First, the maximum available linear extent of a lattice 
is of the order of L = 10 3 , and due to boundary effects, 
only a fraction (let say l/4th of the lattice) can be used. 
As a consequence, the value of r\ a is strongly altered by 
the log-correction in eq. JBJ 3 ) . Keeping a fixed r\ a = \ 
means that we attribute all the deviation of numerical 
data to the existence of Q G , but again, this correction is 
only an effective one, since the log term in the correla- 
tion function J5| is only the first of a series and the next 
term (see e.g. Ref. ^01 Ej) will affect in a similar man- 
ner the value of 6 a and make it an effective one 4 ). A 
second objection comes from the appearence of a typical 
length scale which has strong influence in the correction 
to scale-invariant behaviour. Usually, one does not have 
to take care to the unit length, since the physical quanti- 
ties exhibit scale invariance at criticality. Strictly speak- 
ing, this is no longer true when a correction to scaling 
(e.g. a logarithm) is present, and the above expressions 
for the correlation functions should be rewritten with 
respect to some scale factor a. This factor is related to 
the lattice spacing, but as well to some physics of the 
problem, like the size of vortex pairs in the XY model 
or the typical disorder length in the RBIM, it is thus 
non universal and depends on the model. The local ex- 
ponent becomes 

Ves(r) = rjo- - t-^j- . (7) 
In 77 a 

What is important to notice here is the fact that the 
value of a determines the amplitude of variations of the 
effective exponent n e s(r). According to these observa- 
tions, it is probably already a good result to only pre- 
dict a correct sign for the correction exponent and 
we have no stronger ambition in this paper. 

We propose a different approach which is almost free 
from boundary effects and enable us to work with the 
asymptotic expression of the correlation function in the 
infinite plane. The simulations 5 ^ are performed inside 
a finite system, but the functional expression of the 
correlation function inside such a system is predicted 
by a convenient conformal mapping. This method has 
been applied with success to magnetization profiles in 
the case of the pure XY model ^H] and was extensively 
used in the case of dis ordered Potts models 1201 in two 

3 'With 9cr = ±^ and simulations available up to relative dis- 
tances as large as r = 500, one can at most reach values of r\ a in 
the range 0.23 - 0.27. 

4 ) Using the expression G(r) X r 1 / 4 ~ ln 1/s r X (1 + j^ 1 ^ 1 ) 
given in Ref. 1101 for the BKT transition and again with r ~ 500, 
we get 6Ur) = 6 a + A, — l7, ln ,'" Wi« ~ 0.117. 

° u \ J u 16 m r+(m m r)/16 

5 ) Standard Wolff cluster algorithms. 



dimensions where it was shown to provide quite accu- 
rate results. More problematic is the fact that we apply 
a method which is known to be valid at a really scale- 
invariant fixed point, i.e. in the absence of corrections 
to scaling which break (at the correction level) dilata- 
tion symmetry. In the following, we consider systems of 
reasonable sizes {L up to 256), the asymptotic regime 
r — > oo, G ff (r) ~ r~ r,a , is thus far from being reached 
and the variable lnr only varies within a narrow range. 
Accordingly, eq. J2J might be replaced by an algebraic 
decay with an effective exponent, G a {r) ~ 7-~ I?off , and 
displays scale invariance, at least in the range of dis- 
tances under consideration. 




Fig. 1: a) Conformal mapping of the infinite complex 
plane inside a square w. b) Sketch which shows how 
the boundary conditions (BC) follow from the folding 
("pillow" geometry), c) Example of the profile of the 
correlation function between the upper right corner and 
other points in the square with these particular BC. 

- The upper complex plane 3m ( > is mapped inside 
a square —L/2 < Kew < L, < Smw < L through the 
Schwarz-Christoffel conformal mapping £ = sn ^r^- 

- In order to relate this finite geometry to an original in- 
finite plane (with complex variable z), one may use the 
Schwarz transformation £ = z 1 / 2 which has the effect of 
a folding of opposite edges of the square two by two |?T] 
as shown in figure Q 

There, the effect of the conformal mapping, 

G a {w 1 ,w 2 ) ~ \w'(z 1 )\-^\w'(z 2 )\-^G (r (z 1 ,z 2 ) ) (8) 
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is just to define a rescaled distance variable, called 
k(wi,w 2 ), in terms of which one recovers inside the 
square with these special boundary conditions, a sim- 
ple power law for the correlation function: 

G a (w 1 ,w 2 ) ~ [K(w 1 ,w 2 )}~ Vetf , (9) 
k(w 1 ,w 2 ) = |w / (zi)|^|w'(z 2 )|^|zi - z 2 \, (10) 

where \w'{z)\ = ^ |cn(2Kw/L) dn(2Kw/L) sn(2Kw/L)|" 
and G a {zi,z 2 ) is the correlation function in the plane 
(with z — sn 2 (2Ku>/L)). Here, cm, dnx and snx are 
the Jacobi elliptic functions, L the linear size of the 
lattice, and K ~ 1.58255 is a constant related to the 
aspect ratio of the system. 
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Fig. 2: a) Log-log plot of the correlation function 
G a (w\,W2) of the XY model vs the relative distance 
in the original square geometry for a system of size 32 2 . 
b) Same data plotted now vs the rescaled relative dis- 
tance k(wi, U)2). 



The main advantage of this technique is that one lat- 
tice size L is in principle sufficient (provided it is large 
enough) , since the shape effects are included in the con- 
formal mapping and the method is not much sensitive 
to finite-size effects. The effect of discretization of the 
lattice is only appearent at the scale of a few lattice 
spacings. One more advantage is the fact that all the 
information encoded in the correlation function is used, 
since all the points w inside the square enter the fit (see 
fig. El- Now, as we noticed above, it is necessary to take 
into account the existence of the logarithmic term if we 
want to understand the leading singularity. In order 
to emphasize this comment, let us discuss briefly the 
results presented in figure |3] for L = 32 to 256 in the 
case of XY, 4-state Potts and random-bond Ising mod- 
els. We show the log-log plot of G G {wi,w 2 ) with respect 
to the rescaled distance k(wi,w 2 ) for the three models 



under consideration. One observes the remarkable lin- 
ear regime (on this log-log scale) over the whole range of 
variables, and in particular no boundary effects, as these 
were included in the conformal mapping. Nevertheless, 
while the expected slopes should all be equal to the same 
Tic = i, a deviation from this value is suspected. Then, 
a power law fit leads to leading singularities with expo- 
nents 6 ) r? eff (XY) ~ 0.233(3), r? eff (PM) ~ 0.264(6) and 
7? e ff(RBIM) ~ 0.265(14) (there is a slight variation, de- 
pending on the size) . All these results are in poor agree- 
ment with the exact result 4. This is a clear evidence 
that the logarithmic correction has to be taken into ac- 
count. 
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Fig. 3: Log-log plot of the correlation function inside 
the square, G a (w\,W2) as a function of the convenient 
rescaled variable k(w\,W2) (for sizes from L = 32 to 
L = 256 for the three models under consideration). 

Figure 01 shows a semi-log plot of the rescaled corre- 
lation function /(|zi2|) = Ga-(w\, w 2 ) x [k(wi, w 2 )p" 
against the relative distance \z 2 — z\\ in the infinite 
plane geometry. One observes empirically that the ex- 
pected behaviour /(|zi 2 |) ~ A{\n{\z\ 2 \/ a)) 6 " according 
to eq. (|2J) is in fact linear in this scale, that is 



/(|*i 2 |)~A) + £i In |*12 1 



(11) 



This is coherent with the logarithmic correction pro- 
vided that the dimensionless inverse typical length scale 
a -1 is larger than any of the accesible relative di- 
mensionless distances at the sizes available, a" 1 ^S> 
max\zi 2 \. Under these conditions, the logarithmic cor- 
rection yields 



/(IzHD^OnCa- 1 ))*' ! 



hpj ln|Z12 

6 ) Statistics over the 7 values of size L for each model. 



(12) 
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Tab. 1: Values of the effective exponents deduced from 
fits of the curves shown in figs. [3] and 0] In the case 
of RBIM, the value of Bi/Bo is roughly 3 times larger 
than for Potts (same sign), but strongly fluctuating. 













L 


XY 


Potts 


RBIM 


XY Potts 


64 

128 

256 


0.231 
0.233 
0.234 


0.268 
0.266 
0.260 


0.259 
0.287 
0.255 


0.0087 -0.0075 
0.0077 -0.0081 
0.0072 -0.0074 



i.e. 9 a appears in the ratio B\/Bq = CT /ln(a _1 ). The 
slope B\ in eq. (|llfl is positive if 9 a > and negative 
otherwise. 




Fig. 4: Semi-log plot of the rescaled correlation func- 
tion for XY and Potts models vs relative distance in the 
plane geometry (L = 256) . The slope Bi is positive for 
XY model and negative in the 4-state Potts model. 

As we have mentioned previously, the accessible lattice 
sizes are too small, so that we do not expect any precise 
determination of the correction exponent 9 a . More dra- 
matic is the fact that the 9 a exponent appears in our 
expressions mixed with the non universal length scale a 
from which a reliable value would hardly be extracted. 
Nevertheless, what is shown in Fig. 0] is that a loga- 
rithmic correction is consistent with the data for both 
XY and 4-state Potts models, that the signs of the ex- 
ponents of the log terms are opposite for both models 
and that their absolute values are quite close to each 
other (see table 1). In the case of the disordered Ising 
model, a quantitative analysis is made difficult due to 
the fluctuations introduced by the disorder average, but 
the leading exponent n e f[ was compatible with a nega- 
tive correction exponent 9 ai as expected. 
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